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Our  ultimate  aim  is  to  exploit  the  nearness  to  C-auss ianity  of 
velocity  probability  distributions  in  turbulence,  by  expanding  the  velocity 
field  function  about  the  Gaussian  approximation.  Many  of  the  mathematical 
manifestations  of  our  method  are  so  new,  however,  that  it  is  hard  to  obtain 
physical  insight  into  them.  Hence,  we  have  undertaken  a  preliminary,  pilot 
project  of  a  simplified  nature,  which  usee  instead  of  the  Navier-Stokes  or 
MHD  equations  the  Burgers  one-dimensional  model  equation:1 

du  du  „  d2u  .  . 

3t  +  u  SE  ‘  9  ^2  “  0  < l) 

IHie  expansion  is  in  terms  of  Wlener-Hermlte  functionals 

u(x,t)  a  /  K^X^(x-x^)H^ ^(x^)  dxx+  //  K^(x-x1,x-x2)H^(x1,x..;)dx1  dx2+  •••  (2) 

The  K^,  which  are  implicitly  functions  of  t  =  time,  are  ordinary  (non¬ 
random)  functions.  Integrations  are  over  (-  »,»).  The  ^  are  the 
Wiener-Hermite  functionals  as  defined  in  Reference  k.  They  are  random 
functions,  whose  salient  properties  are  (a)  their  mutual  statistical 
orthogonality,  expressed  by  such  relations  as 

^H(i)  -  0  if  i  +  J,  <(h(1)(xi)  H(1)(x2))>  .  S(Xl-x2),  etc.; 

(b)  is  the  "ideal  random  function"^’-',  or  derivative  of  the  Wiener 
function  ’  ,  whence  the  first  term  of  (2)  is  a  Gaussian  function,  and  the 
later  terms  are  in  some  sense  successively  higher  corrections  to  Gaussianity; 

(c)  the  expansion  is  rigorously  convergent  for  abroad  class  of  random 
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functions,  hence  should  converge  rapidly  whenever  the  velocity  distribution 
is  approximately  Gaussian,  which  is  particularly  the  case  in  the  final 
stage  of  turbulence. 

If  this  expression  for  u(x,t)  is  substituted  into  (i),  one  can 
obtain  from  the  statistical  orthogonality  relations  an  infinite  set  of 
equations  relating  the  K^.  When  the  are  Fourier-transformed 
according  to- the  relation 


C^(j 


-1  (1) 


..£  )  d^- 


•d* 


V 


(3) 


these  equations  become 

K(l)(X)  +V  +  •••  =  0  (4) 

Jt  K(2)(*r*2)  +v(x1+'^)2k(2)(x1,x2)  - 1  (^1+^)k(i)(xl)k(i)(x2)  0  (o) 

|t  +  V  (X1+Jt2+^)2K(3)(X1,X2,  ^ )  - 

^  (WV'k(1)(Vk(^V3)  +  k(i)(*2)k(2)(x3,tcl)  + 
k(i)(jc3)k(2)(x1,^))  +  •••  -  0  (6) 


The  dots  in  equations  (4)  through  (6)  stand  for  an  infinite 
nunfcer  of  terms,  involving  kernels  with  higher  indices;  for  example, 
equation  (4)  contains  a  term  -  2iX(K^-  K^2^),  where  the  inner  product  is 
defined  as 

(K^«  K^)  »  ~  /  K^CX’)k^(%,  X)  dX 

Formalizing  the  assumption  that  the  series  (2)  converges  rapidly, 
we  assign  to  such  inner  products  (subject  to  verification)  an  order  of 


-  3  - 


smallness  equal  to  the  sum  of  the  subscripts.  Thus,  the  terms  given 
explicitly  in  equation  (4)  are  of  first  order,  the  omitted  term  Just  cited 
is  of  third  order,  and  ail  others  are  of  still  higher  order.  If  we  retain 
only  the  terms  written  down  in  equations  (4)  through  (6),  we  have  approxi¬ 
mate  relations,  homogeneous  in  order  of  smallness,  which  determine 
(subject  to  initial  conditions)  each  of  the  first  three  kernels  to  lowest 
order . 

As  initial  conditions  we  asstime 


K(l)(Y)  I.  .  i  AXm  e"1 
1  t«o 


2  2 


k(2V1,*2)  it„  ■  i  B<*1+  v 


_  -i* (*,*+%*) 

i  B »,+  2L)n  e  1  2 


k (3)Ck1,k^,^)  1^  -  i  c (VV  V*  e 


-i2(\2+afe2+^2) 


(8) 

(9) 

do) 


where  A,  B,  and  C  are  constants.  These  conditions  are 'initial"  ones  only 
when  considered  relative  to  a  stage  of  decay  late  enough  to  be  character- 
liable  by  a  single  correlation  length  I.  The  exponentials  are  character¬ 
istic  forms  for  the  late  decay  of  the  dlffuslon-llke  equations  (4)  to 
(6)  .  The  terms  X  ,  (X^+  Xg)  ,  •••  are  presumably  the  leading  terms  of 
Taylor  expansions,  and  a,  n  are  assumed  Integers.  One  finds, that  a  ■  0  leads 
to  (K^:K^)  (double  inner  product)  constant  in  time,  which  means  a 
non-decaying  mean  velocity,  which  is  physically  untenable.  The  choice 
m»  1  leads  to  (K^:IC^2^)  decaying  and  at  the  same  rate  as  (K^*K^). 

As  our  starting  point  we  chose  m  -  n  ■  p  -  1.  We  have  Integrated 
equations  (4)  -  (6)  (we  omit  the  expressions  for  the  K's  so  obtained,  for 
lack  of  space)  and  obtained  some  simply  interpreteble  results  from  them: 
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The  omitted  terms  in  the  equations  of  the  K^,  when  evaluated 

In  terms  of  these  approximate  solutions,  decay,  relatively  to  those  re- 
-l/2 

talned,  as  t  '  or  faster.  This  guarantees  the  self-consistency  of 
omitting  them.  (This  result  is  obtained  by  simple  dimensional  analysis 
of  the  expressions  involved.) 

The  distribution  approaches  a  Gaussian  one  in  time:  the  flatness 
factor  of  the  one-point  distribution  is,  in  terms  of  the  K^,  and  to  the 
present  order  of  approximation, 

T~  3 .  2  fafryg>)-te(i)-*(g)),>  Hi 

(lUJ.jU))* 

The  fraction  in  the  second  term,  apart  from  the  factor  24,  is  a  measure 
of  the  relative  deviation  from  Gauss ianity  ("R.D.G."),  which  we  denote 
by  <f>.  The  present  assumptions  lead  to  <f>  decaying  as  t-1^2.  (It  is  impor¬ 
tant  to  note  that,  contrary  to  what  one  might  think,  this  is  consistent 
with  the  equal  rates  of  decay  mentioned  above  for  (K^:K^)  and  (K^*K^. 

As  a  venture  in  investigating  the  predictive  value  of  the  inherent 
mathematical  structure  of  the  R.D.O.  in  this  scheme,  we  undertook  to  evaluate 
the  R.D.G.  for  space  derivatives  of  the  velocity  field.  As  is  well  known, 
in  turbulence  experiments,  this  quantity  is  positive  and  increases  strik¬ 
ingly  with  the  order  of  the  derivative^.  If  ^  denotes  the  R.D.G.  of  the 
nth  derivative, 


K  •  -  -V 1/2  ♦*4i*  *2  Ko  ♦  » 


24V 


n2  '  “  rny 


where  ^  are  numerical  quantities  obtained  from  the  and 


(12 
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a  *  4  T>  B/A2  -  1 


(13) 


b  »  1  -  12  1>  B/A2  -  2hTp  C/A3 

Computed  values  of  the  J  ■  0,  1,  2,  3  are  given  In  the  table  below: 


n 

**n0 

*n2 

*n3 

0 

4.3 

4,3 

4.5 

-.75 

1 

7.1 

9,5 . 

6.4 

-2.3 

2 

14.3 

8  ,9  .  .  . 

5.9 

. -4,5 

_ 

20.6 

12.8 

6.2 

-10.0 

In  discussing  these  results  we  suppose  that  (l)  has  been  made 

-3/2 

dimens ionless ,  u  being  measured  in  units  of  uq  =  Ai  ,  x  in  units  of  2, 

t  in  units  of  2/uq:  then  T?  is  R**1  where  R  =  ,  the  Reynolds  number 

at  t  »  0.  The  asymptotic  result  given  in  (12)  requires  times  larger  than 

the  Reynolds  number.  Although  it  is  hoped  that  the  representation  used 

here  will  be  useful  in  the  description  of  flows  during  times  when  the 

Reynolds  number  is  large,  we  content  ourselves  here  with  a  regime  where 

the  Reynolds  number  is  small  or  of  order  unity.  In  such  a  case  ~\>  is 

large  or  of  order  unity  and  the  correction  is  small  if  *i2t>?l.  Here 

a  and  b  may  take  on  various  values,  difficult  to  determine  from  the  present 

state  of  the  theory.  In  general,  for  reasonable  choices  of  A,  B,  and  C 

the  deviation  from  Geussianity  increases  with  increasing  order  of  the 

derivative.  For  Instance,  it  is  easily  seen  that  there  is  an  Increase 

o 

if  C  -  0,  and  4  Vb/A  >1.  This  choice  represents  the  situation  where 
at  the  initial  instant  (of  course  after  the  flow  has  decayed  as  described 
above)  ■  0,  and  is  order  unity  or  small  but  not  zero.  This 


(HO 
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increase  in  non-Gaussianity  with  order  of  the  derivative  is  observed 

experimentally  for  real  fluid  flov^  although  for  far  shorter  decay  times. 

It  is  Interesting  to  compare  these  results  with  calculations  which  have 

been  made  of  the  statistical  characteristics  of  solutions  of  (l)  by 

Moomaw^ .  He  finds  that  the  correction  to  the  Gaussian  value  for  the 

-l/2 

flatness  factor  is  proportional  to  t  '  also;  however,  the  dependence  on 
the  Reynolds  nunber  is  different.  In  the  main,  these  results  are  en¬ 
couraging  and  suggest  that  the  stochastic  representation  used  here  may 
be  useful  in  the  treatment  of  turbulence  problems. 
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